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ABSTRACT 

We combine data from the MGC, SDSS and UKIDSS LAS surveys to produce 
ugrizYJHK luminosity functions and densities from within a common, low redshift 
volume (z < 0.1, ^ 71, 000/i^^Mpc'^ for L* systems) with 100 per cent spectroscopic 
completeness. In the optical the fitted Schechter functions arc comparable in shape 
to those previously reported values but with higher normalisations (typically 0, 30, 
20, 15, 5 per cent higher </)*-values in u,g,r,i,z respectively over those reported by 
the SDSS team). We attribute these to differences in the redshift ranges probed, in- 
completeness, and adopted normalisation methods. In the NIR we find significantly 
different Schechter function parameters (mainly in the M* values) to those previously 
reported and attribute this to the improvement in the quality of the imaging data 
over previous studies. This is the first homogeneous measurement of the extragalactic 
luminosity density which fully samples both the optical and near-IR regimes. Unlike 
previous compilations that have noted a discontinuity between the optical and near-IR 
regimes our homogeneous dataset shows a smooth cosmic spectral energy distribution 
(CSED). After correcting for dust attenuation wc compare our CSED to the expected 
values based on recent constraints on the cosmic star-formation history and the initial 
mass function. 

Key words: galaxies: fundamental parameters — galaxies: luminosity function, mass 
function — galaxies: structure — galaxies: statistics — infrared: galaxies — surveys 



1 INTRODUCTION 

The galaxy luminosity distribution at any given epoch is 
a fundamental observable feature of the universe. It is an 
account of how the space density of galaxies varies with 
flux, and provides an insight into how visible matter is 
fragme nted. This can be u sed to constrain both galaxy evo- 
lution |Benson et al.ll200"3 ) and structure formation models 
l|Cole et al.l l200dl . For over half a century astronomers 
have attempted to parametrise the galaxy luminosity 
distribution. The current standard, the Sc hechter f unctio n 
dSchechtej 1 19761 ). built upon work by iHubbld (|l936D . 
IZwickvl (ligSTl ) and lAbeiil (|l958l '). amongst others. 
Observations in visible light, and in particular the wave- 
length range 350 — 550nm, are generally dominated by 
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the most recently formed stellar population, whereas the 
light in the near infrared (NIR) is typically dominated 
by the longer lived, lower mass stars that constitute the 
bulk of the stellar mass (modulo some contamination from 
the AGE branch). The N IR is also less susc eptible to 
internal dust attenuation (|Calzetti et al.l 1 19941 ). with an 
estimated ~ 80 per cent, ~ 50 per cent and ~ 20 per cent 
of the integrated flux from galaxies being a ttenuated in 
the UV, optical and NIR respectively (see [Driver et al.l 
|2008| ). These two benefits, the focus on stellar mass and 
lower attenuation, make NIR luminosity functions arguably 
more useful for characterising the underlying properties of 
the galaxy population, for comparison with semi-analytical 
models, and, as a stepping stone to calculating the yet 
more fundamental stellar mass function (see for example 
iBaldrv et al.|[200^ ). However recent results, particularly in 
the K band, show a relatively large range in the reported 
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Schechter function parameters and do not yet span the full 
NIR wavelength range available; for example no Y band 
LF has yet been reported. The first aim of this paper is 
to provide fresh estimates of the NIR luminosity functions 
using the latest available data from the UKIRT Infrared 
Deep Sky Survey Large Area Survey in the Y, J,H & K 
passbands. 

By integrating over the product of luminosity and the lu- 
minosity function, one derives another useful measurement: 
the total luminosity density, j (sometimes written as £) 
at a specific wavelength. The form this takes for a single 
component Schechter function in solar luminosity units 
(/iLq ;^Mpc~^) can be seen in Equation [1] 



(1) 



and a\ are the wavelength dependent 
Schechter function parameters. Measuring luminosity den- 
sities across the full UV/optical/NIR wavelength range 
(where starlight entirely dominates the energy output), al- 
lows one to build up the cosmic spectral energy distribu- 
tion (CSED) for a representative volume of the local Uni- 
verse. This is important as it provides a description of the 
mean radiation field from UV to NIR in which the con- 
tents of the nearby Universe are lying. However the ra- 
diation we detect is only some fraction of that produced 
with a signific ant amount atten uated by dust within the 
host galaxies (iDriver et al.ll2007^. C orrecting for the dust 
attenuation (see Driver et al.l 20081 ) . one can derive the 
pre-attenuated CSED. We can use this empirical result to 
test our underst anding of th e cosmic star-formation his- 
torv (e.g. iMadau et al.l Il998l , iBaldrv fc Glazebrookl |2003| . 



IWilkinsetal .1120081 ) 
One of the reasons this approach has not taken off is an 
apparent discontinuity between the optical and NIR lu- 
minosity density measurements which i s clearly unphysi- 
cal. This was first reported in IWrighd i 200 ih where an 
extrapolation of the SDSS Early Data Release luminos- 
ity density measurements in ugriz led to an overpredic- 
tion of the observed NIR luminosity densities by a fac- 
tor of 2-3. This issue was mainly resolved when revised 
SDSS LF estimat es were published wh ich incorporated evo- 
lutionary effects (jBlanton et al.l l2003l ) , however the prob- 
lem has not entirely v anishe d and is clearly noticeable in 
iBaldrv fc Glazebrookl (|2003h where data from t he SDSS 
( Blanton et al.ll2003l). GALEX (IWvder et al.ll2005l) a r id the 



NIR studies of ICole et al.l (|200ll ). iKochanek et all (|200ll ') 
were used to produce the nearby CSED. This analysis actu- 
ally brought to light a number of issues including the incom- 
plete coverage of the NIR wavelength range (with only J and 
K values available), the significant scatter in the available 
-ff-band data, and an apparently unphysical step-function 
between the optical and the J-band resulting in the J-band 
data being dropped. This di scontinu i ty fro m the optical to 
near-IR echoes that seen in IWrighd i 200 ll ) and is also ob- 
vious when comparing the recent 6dfGS surv ey measure- 
ment s of the galaxy luminosity density in JHK (* Jones ct al. 
200(: ) with recent SDSS results i n ugriz (jBlanton c t al. 



20031 : iMontero-Dorta fc Pradall2009l ) . In general any discrep- 



face brightness selection bias in the NIR data (|Smith et al.l 
l2009l) . discrepancies in the photometric calculation, or spec- 
troscopic incompleteness bias. Certainly, when one reviews 
the most recently published NIR luminosi ty densities one 
does fi nd s ignificant scatter. I n particular, iKochanek et al.l 
(|200ll ) and lHuang et al.] (|2003l ) examined insufficiently large 
vo lumes to ov e rcome cos mic variance. Other attempts, such 
as ICole et all (|200ll ) and Ijones et all (|2006l '). have probed 
greater volumes but are dependent on shallow 2MASS imag- 
ing data that has been shown to be susceptible to sur- 
face brightness (SB ) bias, missing b o th galaxies and flux 
fe.g. lAndreonI l2002l . iBeU et al.l |2003| . lEke et al. I l2005l and 
iKirbv et al.ll2008h . 

The UKIRT Infrared Deep Sky Survey Large Area Sur- 
vey (UKIDSS LAS) extends ~ 2.7 mag deg"^ deeper than 
2MASS m K, and the r efore should be less susceptible to 
SB effects. ISmith et all (|2009l ') recently produced a if-band 
SDSS-UKIDSS luminosit y function that appe ars to agree 
closely with the results of lCole et al.l (|200ll ) and I Jones et al.l 
(|2006l ). This would suggest that 2MASS and UKIDSS pho- 
tometry are consistent and the reported NIR LFs robust. 
However, due to an unresolved issue with the UKIDSS ex- 
traction software affecting ^ 10 per cent of the UKIDSS 
data (see detailed discussion in Appendix El). ISmith et al] 
question the validity of their own results. ISmith et al. also 
restricted their analysis to K-h&xiA only where, for the pur- 
pose of recovering the cosmic SED, it is obviously desir- 
able to recover measurements in all available filters (i.e. 
Y JHK). Fin ally a significant issue may well arise from cos- 
mic variance (|Driver fc RobothamlboiOl ') with NIR surveys 
typically being based on more local samples (e.g., 6dfGS 
with {z) ~ 0.05) and the optical data being based on deeper 
samples (e.g., SDSS with {z) ~ 0.1). In an ideal situation 
one would seek to derive the full cosmic SED from within 
a single survey where the impact of cosmic variance will af- 
fect the measurements at all wavelengths in a similar man- 
ner (modulo colour dependent clustering). In this paper we 
do precisely this; by combining data from the Millennium 
Galaxy Catalogue (MGC), SDSS and UKIDSS LAS surveys 
we derive the ugrizY J H K luminosity distributions and pre- 
and post attenuated cosmic SEDs within a single well un- 
derstood volume of ~ 71, 000/i^Mpc^. 

In Section [2] we briefly describe the three surveys, their 
shared area and appropriate flux limits (to ensure against 
colour bias). In Section |3] we describe our methodology 
for deriving Schechter luminosity function parameters. Fi- 
nally, in Sections U and [S] we present our results, discuss 
how they compare with previous work and what conclu- 
sions we can make. We also include an outline ( Appendix 
of the problems with the CASU source extraction in early 
UKIDSS data releases, and the steps we undertook to re- 
analyse UKIDSS data for this paper. Throughout we adopt 
axih = l,Ho = lOO/i km s~^ Mpc"\ Q.m = 0.3, = 0.7 
cosmological model. All optical magnitudes are quoted in the 
SDSS photometric system (which is consistent with the AB 
system ±0.04mag) and all NIR magnitudes in the UKIDSS 
preferred Vega system unless otherwise stated (Table. 1 
shows the relevant conversions between the AB and Vega 
systems) . 



ancy in the CSED between optical and NIR wavelengths 
could in part be explained by: a very top heavy IMF, c osmic 
variance in the NIR data (e.g. ISomerville et al.]|2004l ). sur- 
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Figure 1. The distribution of galaxies along the common MGC- 
SDSS-UKIDSS region for (top-to-bottom) the uBgrizYJHK 
bands (as indicated). Note that the MGC data does not pro- 
vide contiguous coverage within the designated area due to gaps 
between CCDs (the "L" shaped layout of the CCDs causes this 
problem) and masked regions around bright stars and satellite 
trails. The SDSS DR5 coverage of the MGC is complete while 
the UKIDSS LAS shows some clear gaps. The coverage as a frac- 
tion of the total area is summarised in Column 2 of Table [2] 

2 DATA 

In this paper we combine data from three distinct datasets : 
the Millennitim Ga laxy Catalogue (MGC: iLiske et al.|[20o3 : 
Driver et all 1200511 . the Sloan DiRital Sky Surve y (SDSS; 
York et al.ll200Gl : lAdelman-McCarthv et al.ll2007l ). and the 
UK I nfrared Deep Sky Survey Large Area Surv ey (UKIDSS 
LAS: iLawrence et akHioOTl : IWarren et al.ll2007[ ). The MGC 
provides the deepest {B < 20 mag) wide area redshift survey 
to date (extending about L5 mag fainter than SDSS main), 
while the SDSS provides full optical coverage (ugriz), and 
the UKIDSS LAS fuU near-IR coverage (YJHK). The com- 
bination of these datasets enable the full NIR-optical map- 
ping of the cosmic energy output of the Universe emanating 
from a representative volume, robust to both wavelength de- 
pendent cosmic variance and spectroscopic incompleteness. 
In this section we briefly introduce each survey and describe 
the coverage and flux limits of the common region. 

2.1 The Millennium Galaxy Catalogue (MGC) 

The Millennium Galaxy Catalogue (|Liske et al.ll2003l ) is a 
deep {Biim = 24 mag, //lim — 26 mag arcsec ), B band 
galaxy survey created using the Wide Field Camera on the 
2.5m Isaac Newton Telescope, with observations taken be- 
tween 1999 and 2001. Surveying a long (75°), thin (0.5°), 
equatorial strip amounting to 37 deg'^ of sky (30.88 deg^ 
after cleaning and cropping), the MGC contains data on 
over a million objects, with 10,095 galaxies brighter than 
B = 20 mag (this resolved sub-catalogue is referred to 
as MGC-Bright, and the integrity of every object within 
it has been verifled by eye and fixed where necessary). 



Using the extensive MGC overlap regions, for objects in 
the range 17 ^ -Bmgc ^ 21 mag, the astrometric rms 
uncertainty has been shown to be ±0.08 arcsec, and the 
internal photometric uncertainty ±0.023 mag. The seeing 
ranged from 0.9 to 2 arcsec, with a median of 1.3 arcsec. 
Object det ection was achieved usi ng the SExtractor pro- 
gram (see iBertin fc ArnoutsI Il996h . with a constant sur- 
face brightness threshold, /inm, of 26 mag arcsec"^, and 
galactic extinc t ion c or rected using dust maps created by 
ISchlegel et all (|l998l ). [Driver et all l|2005l ). have obtained 
redshifts for 96 per cent of the MGC-Bright galaxies us- 
ing the SDSS and 2dFGRS d ata releases combi ned with a 
dedicated MGCz 2dF survey. lAUen et al] (|2006l ) have con- 
ducted bulge disc decomposition using GIM2D for the entire 
MGC-Bright catalogue. All data are publicly available from 
htt p : / / www, eso . org/ ~ j liske /mgc /] 

2.2 The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS. lYork et al.ll2000l ) is 
the largest combined photometric and spectroscopic survey 
ever undertaken, containing spectra of ^930k galaxies and 
imaging of over 11500 deg'^ of sky, using five filters with av- 
erage wavelengths between 300 and 1000 nm (ugriz). SDSS 
data has been publicly released in a series of 7 data releases. 
We are using the fifth data release catalogue (SDSS-DR5), 
which covers approximately 8000 deg'^ and contains spec- 
troscopy for 675k galaxies. The MGC region falls within the 
SDSS DR5 area of coverage. 

Using the SDSS DR5 database stored in the WFCAM Sci- 
ence Archiv^B, we output a complete list of SDSS Photo- 
Primary objects within the MGC region. We then matched 
MGC-Bright objects to SDSS Photo Primary objec ts using 
the STILTS catalogue matching tool (|Tavlodl2006l ). with a 
maximum centroid separation tolerance of 2.5 arcsec. Where 
there are multiple matches within 2.5 arcsec to one MGC ob- 
ject, STILTS takes the closest matching object. For the pur- 
poses of photometry we adopt the extinction corrected SDSS 
Petrosian apparent magnitudes (PetroMag_y — extinction^ , 
where X is u, g, r, i or z). Our final DR5-MGC-Bright 
matched catalogue contains 10050 SDSS-MGC matching 
galaxies (99.6 per cent of the 10095 sources that make up the 
MGC-Bright catalogue). The coverage of the MGC by SDSS 
DR5 is shown in Fig. [T] SDSS fiuxes are reported through- 
out in AB magnitudes, conversions to Vega are shown in 
TablelU 

2.3 UKIRT Infrared Deep Sky Survey (UKIDSS) 

UKIDSS (|Lawrence et al.ll2007l ) is a seven year near-infrared 
survey programme that will cover several thousand deg'^ of 
sky. The programme utilises the Wide Field Camera (WF- 
CAM) on the 3.8m United Kingdom Infra-Red Telescope 
(UKIRT). The fuU UKIDSS program consists of five sep- 
arate surveys, each probing to a different depth and for a 
different scientific purpose. The shallowest of these surveys, 

1 The WSA l lHamblv et al.ll2008l) is a storage facility that con- 
tains copies of catalogues from a number of surveys, including the 
SDSS and the MGC, as well as being the primary store for raw 
UKIDSS data. 
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Table 1. The Absolute magnitude of the Sun in various fil- 
ters for the AB and Vega systems along with the approximate 
filter central wavelength. These values were derived for us by 
Paul Hewett (priv . comm ). Mq differ from those in Table 1 of 
iBlanton fc Roweij bOO?! ') by mag, 0.03 mag, 0.07 mag, 0.03 mag 
and 0.03 mag in the u, g, r, i and z bands. 

TJpm) Mq{AB) MQ(Vega) Mq(AB) - MQ{Vega) 



u 


0.3546 


6.38 


5.47 


0.91 


g 


0.4670 


5.15 


5.23 


-0.08 


r 


0.6156 


4.71 


4.55 


0.16 


i 


0.7471 


4.56 


4.19 


0.37 


z 


0.8918 


4.54 


4.00 


0.54 


Y 


1.0305 


4.52 


3.89 


0.63 


J 


1.2483 


4.57 


3.63 


0.94 


H 


1.6313 


4.71 


3.33 


1.38 


K 


2.2010 


5.19 


3.29 


1.90 



SDSS DR5 provides complete coverage of the cleaned 
MGC region in all bands (see Fig. [T] ugriz). A detailed 
comparison of the photometry, astrometry, and deblending 
between th e MGC and DRl was described in detail in 
ICross et al.1 (|2004l '). In total MGC-Bright contains 10,095 
galaxies to Bmgc = 20.0 mag of which matching detections 
are identified in the SDSS survey for 99.6 per cent of the 
MGC targets, the majority of failed matches comes from 
extreme low surface brightness systems and differences in 
deblending decisions. 

The UKIDSS LAS DR3PLUS does not have complete 
coverage of the combined MGC-SDSS region (Fig. [TJ 
YJHK) and the level of completeness varies for each 
filter. The gaps are due to data failing the UKIDSS data 
control process (typically seeing and sensitivity critereon). 
From the common regions the distribution of galaxies that 
have been flagged as containing major errors or requiring 
deblending (see Appendix [XJ, show no obvious bias along 
the strip (as one would expect given the quality control 
process). The area of overlap can therefore be derived 
based on the fraction of MGC galaxies with UKIDSS data 
available. The coverage of the SDSS and UKIDSS-LAS 
datasets is summarised in Columns 2 and 3 of Table [2] 
The MGC is taken as the master catalogue in all that 
follows as it is the deepest (in terms of flux sensitivity), 
contains the highest signal-to-noise detections and has been 
fully masked and eyeballed. Also, where necessary, all ob- 
jects within it have already been reconstructed or deblended. 



2.5 ugrizY J H K magnitude limits 

In order to derive luminosity distributions for multi- 
wavelength data from a B band spectroscopic sample it is 
important to consider the colour bias. Figure [2] shows the 
Bmgc — X y X colour plots where X denotes ugrizYJH or 
K. The most conservative approach to deriving an unbiased 
luminosity distribution is to simply define a complete sam- 
ple within each band, i.e. cut the sample at a sufficiently 
bright flux where the full colour distribution is fully sam- 
pled (see Fig. (2] long dashed line) and where the number- 
counts have yet to show any indication of a turn-down 
(Fig. 131). Fig. [3] also shows that we can extend our s ample s 
deeper than those used in iMontero-Dorta fc Pradal Hooi); 
our turn-down occurs roughly one magnitude deeper in each 
SDSS passband. We show the iMontero-Dorta fc Pradal lim- 
its on Fig. [2] as long dash-short dash lines. In all SDSS pass- 
bands they are more conservatively cut than our most con- 
servative limit, particularly in the u band. The most lib- 
eral approach is to use all the data and deflne a unique 
flux limit for each individual object based upon the spec- 
troscopic limit of Bmgc = 20.0 mag combined with the ob- 
ject's colour, i.e. Xumit = 20.0 — {B — X), using these B 
and X limits to appropriately weight each object (Fig. [5J 
dotted line). While the former will reduce the sample size 
significantly the latter will incorporate a large quantity of 
data in the regime where flux measurements may not be 
credible. Here we adopt a hybrid approach (Fig. [2] short 
dashed line) where we determine the mean colour for each 
sample (solid line) and combine this with the spectroscopic 
limit of Bmgc = 20.0 mag to determine a nominal limit, 
i.e. the limit where the mid-point of the colour distribution 



the UKIDSS Large Area Survey (LAS), contains the full 
MGC region. 

When complete, the LAS will cover 4000 deg^ of sky to tar- 
get depths (5(7 point source detections in Vega) of K = 18.2 
mag, H — 18.6 mag, J = 19.9 mag (after two passes; this 
paper uses only the first J pass which is complete to 19.5 
mag) and Y = 20.3 mag (for conversion to the AB system 
pleeise see the offsets shown in Table. (TJ. UKIDSS data are 
required to have a seeing FWHM of < 1.2 arcsec, photomet- 
ric rms uncertainty of < 0.02 mag and astrometric rms of 
< 0.1 arcsec. Each position on the sky is targeted for 40s 
per pass. All survey data for this paper is taken from the 
third data release (DR3PLUS). 

During the course of our analysis a number of problems 
were encountered with the online catalogues generated by 
the Cambridge Astronomical Survey Unit (CASU), partic- 
ularly the calculation of Petrosian fiuxes fo r deblended sys- 
tems (as first noted in ISmith et al.l l2009l l . The problems 
are described in detail in Appendix |X] and were considered 
sufficiently insurmountable to warrant re-deriving the cat- 
alogues in entirety from the reduced data. The derivation 
of our final MGC-Y,J,H and K catalogues using SExtractor 
l|Bertin fc Arnoutd[l996l ') are also described in Appendix 1X1 
In summary after revision of the UKIDSS photometry we 
estimate that individual fiuxes in the Y, J, H, and K bands 
are credible to ±0.05 mag. This value derives from a com- 
parison between the original UKIDSS LAS data for which 
deblending was not required, and our revised catalogue. 

2.4 The MGC-SDSS-UKIDSS LAS common area 

Figure IH shows the coverage of the MGC region by the 
SDSS DR5 and UKIDSS LAS datasets. Note that the MGC 
footprint itself does not provide continuous coverage within 
a rectangle but rather a square-tooth profile (the INT WFC 
has a thick "L" -shaped footprint). The MGC region has 
also been carefully masked to remove objects close to bright 
stars, where the fiux might be compromised, satellite trails, 
CCD defects, and edge effects within the MGC region 
which reduces the effective area on sky from ~38 deg^ to 
30.88 deg^ This process negates any issue of the blocking 
factor of distant galaxies by bright foreground stars as well 
as false detections due to spurious light scattering (e.g. 
diffraction spikes and ghosting). Fu ll details of the M GC 
footprint and its masking is given in iLiske et al.1 (|2003l ). 
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Figure 2. The colour-magnitude diagrams for the ugrizY J H K filters versus B. The vertical lines are our apparent magnitude limits in 
ugrizY J H or K, the diagonal linos are our B band magnitude limits. 



Table 2. Parameters defining the coverage and depth(s) of the joint MGC-SDSS-UKIDSS common region along with the adopted K(z) 
corrections and E{z) ranges. Though we set E{z) = 0, we use the /3 ranges in column 8 to calculate the scale of the uncertainty due to 
evolution, and the range of _R'(2)-corrections in column 7 to calculate the uncertainty due to the ftT-correction. The sample size column 
gives the number of galaxies brighter than Limit 2 within the defined redshift limits. The Am^—o i column gives the range of effect the 
combined K+E correction can have on a z = 0.1 galaxy. 



Filter 


Coverage 
(per cent) 


Area 
(deg2) 


Limit 1 
(mag) 


Limit 2 
(mag) 


Sample size 
(0.0033<2<0.1) 


k{z) 


E{z) f} range 




=0.1 


u 


100% 


30.88 


20.0 


20.84 


3267 


(2.16^15:^7)^ 


-1.36—0 


-0.01- 


-0.31 


9 


100% 


30.88 


19.2 


19.61 


3328 




-0.68—0 


0.09- 


-0.35 


r 


100% 


30.88 


17.9 


18.76 


2781 


(0.95l°J*)z 


-0.45-0 


0.00- 


-0.13 


i 


100% 


30.88 


17.3 


18.34 


2623 


(0.48l°-J|)z 


-0.34-0 


-0.01- 


-0.09 


z 


100% 


30.88 


17.0 


18.07 


2437 


(0.03l°;?|)2 


-0.27-0 


-0.01- 


-0.04 


Y 


77% 


23.69 


16.3 


17.38 


1798 


(0.00l°;?*)z 


-0.23-0 


-0.04- 


-0.02 


J 


81% 


25.07 


16.0 


16.89 


1589 


(-o.6iiS:?o)- 


-0.19-0 


-0.08— 


-0.03 


H 


91% 


27.89 


15.0 


16.12 


1890 


(-0.28t«:?^)z 


-0.17-0 


-0.06- 


-0.00 


K 


91% 


27.99 


14.5 


15.67 


1785 




-0.15-0 


-0.16— 


-0.13 
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Figure 3. The number of galaxies, by apparent-magnitude in the 
ugriz Y J H K filters. The dash - dot-da s h lines are galaxy counts 
from iMontero-Dorta fc Pradal | |2009|) . iMontero-Dorta fc Pradal 
justify their early downturn as an issue of redshift incomplete- 
ness, and introduce conservative magnitude cuts accordingly. 



Table 3. Median colours and Scr clipped standard deviations 
above the completeness limits defined by Flux limit 1 in Table. [21 



Colour 


Median 


Std (3(T clipped) 


u- B 


0.81 


0.40 


B~9 


0.38 


0.13 


B - r 


1.21 


0.44 


B - i 


1.62 


0.57 


B - z 


1.94 


0.59 


B-Y 


2.56 


0.62 


B~J 


3.06 


0.51 


B - H 


3.82 


0.59 


B ~ K 


4.25 


0.59 



is spectroscopically sampled. Each galaxy is then allocated 
a flux limit which is the brighter of the nominal limit or 
that defined by the locus Xumit = 20.0 — (-Bmgc — X), see 
short dashed lines on Fig. [2] The conservative limits (Limit 
1), nominal limits (Limit 2), and effective sample sizes are 
shown in Table [S] Table [3] shows the median colour and 3a- 
clipped standard deviation derived from the data above the 
conservative limit (Limit 1). After imposing our magnitude 
limits, every galaxy within our samples has a redshift. Fig- 
ure [4] shows the distribution of our apparent magnitude cut 
samples by redshift. 



— = m) = 0.4 mio r ^ ^,oo..(M--i, (2) 



Where M* is the characteristic magnitude above which the 
power law behaviour is suppressed, with the power law hav- 
ing a slope of Of + 1. The </>* parameter is a normalisation 
constant. Whilst only an approximate match to the galaxy 
luminosity distributior(3, it is the standard parametric fit 
adopted in the literature and therefore useful for comparing 
with earlier work. 



3.1.2 SWML with variable magnitude limits 

There are a number of different techniques for measuring 
the galaxy lum inosity distribution (see the comparison by 
IWillmei|[l997l '). We adop t the SWML m e thod, presented 
and described in detail in lEfstathiou et al] (|l988l ). This is a 
maximal-likelihood method of calculating volume corrected 
binned luminosity distributions in a non-parametric man- 
ner. The resulting luminosity distribution can then be fit 
via x'^-minimisation to determine the Schechter function pa- 
rameters. As noted, the SWML method makes no a priori 
assumption of the form of the luminosity function, but does 
require an additional method for normalisation. Here this 
is done by (1) selecting an absolute magnitude range (Mi 
to M2), (2) finding the number of galaxies inside this range 
{N\^2), (3) calculating the maximum and minimum red- 
shifts over which galaxies of magnitude Mi to M2 can be 
seen and (4) integrating the standard expression for the 
volume interval (Vi^2 = ^) over this redshift range for 
the specified cosmology. The SWML luminosity distribu- 
tions are then scaled such that: </i(Mi to M2) = -^J^^- 
The normalisation is therefore a 1/Vmax-method with the 
SWML luminosity distribution rescaled to produce the re- 
quired number of galaxies within the specified absolute mag- 
nitude range. Care must be taken to ensure that the cali- 
bration volume i s comple t e for the range of absolute mag- 
nitudes selected. IWillmeil (|l997l ') has found that, in samples 
where the faint end of the luminosity distribution is under- 
represented, the behaviour of the SWML method can be 
eccentric, with a larger a than expected being recovered. In 
lEfstathiou et al.] (|l988f ) the method adopts a constant flux 
limit for all galaxies. We introduce a minor modification to 
accommodate for the non-uniform flux limit (because of the 



3 METHODS 

3.1 Luminosity distribution and function 
measurement 

3.1.1 The Schechter function 

The commonly adopted functional form of the luminos- 
ity function is the Schechter luminosity function (jSchechteil 



^ I Jones et al.l ||2006|'| found th a t it d oes not turn down sharply 
enough at M* . iBlanton et al. I ll2003l'l found that there is a de- 
viation from the Schechter function at the luminous end and a 
stro ng positive correl ation between the M* and a parameters, 
and iBell et ahl ||2003| ) have attempted to modify it to fit their 
dataset more accurately. 
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colour bias) with our final statistic given in Equation [S] as: 



<f>{Mk)dM 



W{x) 



where: 

dM 
2 

otherwise 



if - ^ ^ :r s; ^ 



(3) 



if a; < 



dM 
2 



dM 
2 

otherwise 



x_ - f _ dM_ ^ ^ ^ dM_ 

dM " 2 ^ ^ 2 



Tlie key distinction here is the introduction of an individual 
magnitude limit (miim.i) for each object. This modification is 
identical to that adopted by the 2dFGRS, 6dfGS and MGC 
teams when dealing with non-uniform flux limits. 




3.1.3 Calculation of absolute magnitudes 

As the SWML method in Equation [3] requires absolute mag- 
nitudes it is first necessary to convert from apparent to ab- 
solute magnitudes. This requires a cosmological model and 
K{z) and i5(i;)-corrections. The adopted cosmological pa- 
rameters are Oa/^O- S, Q,k=Q.7, Ho=1Q Q k m s~^ Mpc~^, 
compa tible with both ljones^ et all (|2006l ) and lBlanton et al.l 
l|2003h . We reject galaxies with redshifts outside the interval 
0.0033 ^ z ^ 0.10. The lower limit is defi ned as that require d 
to overcome the local velocity field (see [Driver et allboOSh . 
the upper limit is chosen to minimise the uncertainty inher- 
ent in the adopted K{z) and E{z) corrections and ensure 
a uniform survey volume across all wavelengths. For both 
corrections we elect to adopt global rather than individual 
values and Monte-Carlo over a suitably broad range of un- 
certainty to ensure the final uncertainties on the Schechter 
functi on p arameters is realistic. Un l ike th e iBlanton et al.l 
(|2003l ) and iMontero-Dorta fc Pradal (|2009l ') papers, we cal- 
culate our absolute magnitudes with the SDSS passbands at 
z=0, rather than at z=0.1. 

We calculate the A'(2;)-corrections required for our data (col- 
umn 7 of Table [2| via the 13.2Gyr Sa-type galaxy spectra 
from the synthetic library of lPoggiantil (| 19971 ) , and the range 
of possible i4'(2)-corrections from the 7.2 Gyr Sa-type (the 
lower limit of column 7) and 15.0 Gyr El-type galaxy spec- 
tra (the upper limit). 

We do not believe evolution occurs within 0. 0033 < z < 0.1 
(the first four black points in Figure 11 of iPrescott et all 
I2OO9I show no evidence of any occurring); we note that its ef- 
fects would predominantly impact the M* parameter, with 
only a small change in the a and (jy* parameters. For in- 
stance, using the evolution in the u band luminosity density 
presented in iPrescott et al.l (|2009l ) (/3 = —1.36 using our 
sign c onvention) and the prescription of lPhillipps fc Driveij 
(| 19951 ) given by Equation [^: 



i5(2) = 2.5/31ogi„(l + z) 



(4) 



where z is the redshift of the corrected galaxy, the inclusion 
of an evolutionary correction modifies the best-fitting M* 

^ The /3 par a meter in Equation |4] and the (3 variable in 
IPrescott etal] 1I2OO9I ') will, by definition, have opposing signs. 



Figure 4. The redshift distribution of galaxies within our 
apparent-magnitude limited matehed samples, in bins of z=0.006. 
The dashed lines signify our sample redshift limits. 



parameter by 0.09 mag, the best-fitting a parameter by 0.02 
and the best-fitting 0* parameter by 0.0003 h^ Mpc"''. Our 
quoted results do not employ an i5(z)-correction; however 
we calculate the effect that including one would produce. 
We calculate the uncertainty due to evolution by deriving 
the best-fitting Schechter parameters for five equally spaced 
P values from the range shown in column 8 of Table [2j 
and take the standard deviation in these results as our 
evolutionary uncertainty. The effects of the K{z) -\- E{z) 
corrections are, unsurprisingly, strongest in UV (±0.05 
mag in M* in the u band), an d limited in the NIR (±0.01 
mag in M* in the K band). iBlanton et al] (|2003l ') have 
used the evolution correction as a reason for a flattening of 
the faint-end of the luminosity function; within our small 
redshift range we do not believe that this is the case. 
Finally, luminosity uncertainty can move a galaxy into the 
wrong absolute-magnitude bin, and this effect increases the 
counts within the brighter bins by a greater fraction than 
the fainter bins; i.e. a classical Malmquist bias. Where there 
is a large uncertainty in luminosity this must be compen- 
sated for. However, the typical luminosity uncertainty in 
our data (0.04 mag in K - see section lAll 0.05 mag in H, 
0.05 mag in J and 0.04 mag in Y) is small enough to make 
this unnecessary. 



3.1.4 Normalisation 

Within 30.88 deg^ of sky, we calculate the normalisation 
volume over the redshift range (0.023 < z < 0.097) to be 
71069/i"^Mpc^ The same redshift range is adopted for all 
nine filters to insure against cosmic variance. The redshift 
range is selected to be complete for a 1 mag range around the 
M* point for each of the nine filters. We calculate the scale 
of the global overdensity using SDSS-DR7 spectoscopic cat- 
alogue data for a 5150 deg^ rectangular region of sky (130 to 
236 deg RA, to 58 deg Dec, within our redshift range). We 
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find an average source density of 13.00 de-reddened, M*±0.5 
mag galaxies per deg'^ within this region. In the MGC area, 
we find 15.90 de-reddened, M* ±0.5 mag galaxies per deg'^. 
We assume that the area lost due to bright star holes within 
the SDSS region is the same fraction as that lost in the 
MGC region (~3%). When this is taken into consideration 
the MGC is 19.3% overdense. This overdensity is visible in 
Figure [3] For each filter we modify this volume to take into 
account the variation in the area of coverage, spectroscopic 
incompleteness (for all of our samples this — 1), and the 
global over-density of the MGC region, i.e.: 



Vx = Vs./a,x./c,x//mgc 



(5) 



Where X is the filter, /a.x is the fraction of the MGC cov- 
ered by filter X, fc,x is the sample completeness in filter 
X, and /mgc is the global over-density of the MGC. Using 
these corrected volumes, and the number of galaxies within 
a one magnitude range that contains the M* mag galax- 
ies we calculate the galaxy number-density and scale the 
unnormalised number-densities ((^-values) derived from the 
SWML method to reproduce this object density: 

(t,{M) = NxiM* - 0.5 < M < Af* + 0.5) /Vx (6) 



4 RESULTS AND DISCUSSION 

4.1 ugrizY J H K luminosity distributions, functions 
and densities 

The resulting luminosity distributions and fitted Schechter 
functions for all samples (ugrizYJHK) are shown in Fig- 
ure [5] and the Schechter parameters tabulated in Table |4] In 
general the luminosity functions are reasonable fits to the lu- 
minosity distributions based on the reduced-x'^ values. The 
most notable exceptions are in the z-band, which appears to 
show a tentative upturn at fainter magnitudes (M^ — 17 
mag), and the i-band, which appears to shows the same ef- 
fect fainter than Mi ~ —16.5 mag. We have checked that 
the faint-end upturn is not dependent on the faint end mag- 
nitude limit; it remains if we cut our samples using our con- 
servative limits (Limit 1 in Table [2]), or if we use the much 
brighter limits of the SDSS samples. However, we note that 
the faint-end upturn is confined to these two passbands and 
does not seem to be a general characteristic of all our lu- 
minosity functions. There is no sign of any obvious excess 
of very bri ght systems (possibly du e to the defined redshift 
interval as iMontero-Dorta fc Pradal found an overdensity at 
z > 0.1) and in general the Schechter function provides good 
fits around the knee and brighter. Overlaid on the diagram 
are selected recent measurements from other groups. In al- 
most all cases our results lie above those reported from the 
much larger SDSS survey results. As we are using the SDSS 
photometry, this cannot be a photometric issue. Further- 
more as we have compensated for the MGC's over-density 
by calibrating to the full SDSS it is also cannot be due to 
cosmic variance across the sky. A further possibility is simply 
the difference in adopted K{z) and E{z) corrections as it is 
noticeable that the vertical offset does appear to show some 
wavelength dependence with the offset maximum in g (~30 
per cent) and dropping to a minimum in z (~5 per cent). 
However we note that the MGC data extends approximately 
1 mag deeper in all filters (the SDSS sample limits reported 



in iMontero-Dorta fc Pradal are 19.00, 17.91, 17.77, 17.24, 
16 .97 mag in ugr iz respectively, q.v. Table 2, column 5, and 
in iBlanton et al.l the ugriz limits are 18.36, 17.69, 17.79, 
16.91 and 16.50 mag, q.v. Table 1, column 2), our samples 
have 100 per cent redshift completeness, and we use an iden- 
tical redshift ra nge for all filters (0.023 < z < 0.097). The 
recent study bv lMontero-Dorta fc Pradal (|2009l ) by compar- 
ison has a median redshift completeness of 85 per cent (re- 
ported to be both wavelength and flux dependent) that may 
be partly due to the ~55 arcsec minimum fibre proximity 
of the SDSS spectral survey and, although having signifi- 
cantly brighter fiux limits, is used to probe to significantly 
higher red shifts (z ^ 0.2). Moreo ver, t he normalisation 
adopt ed in iMontero-Dorta fc Pradal is the iDavis fc Huchral 
(|l982t ) method. This uses the entire data set and tends to 
overly weight the normalisation towards the higher redshift 
range where incompleteness may be most severe. Without 
reanalysing the SDSS data using our methodology it is not 
possible to ascertain the exact cause of the normalisation 
discrepancy but it is plausible that it is related to this nor- 
malisation method and any bias in the redshift incomplete- 
ness. 

In the NIR the discrepancies are more dramatic, perhaps as 
expected given the rapid development of NIR technologies. 
We are not aware of any published Y band luminosity func- 
tions for comparison. Once again the offset tends to be that 
our results produce a higher space-density of galaxies. This 
is perhaps expected given the significa ntly deeper i maging 
data. Comparing our K-hand result to ISmith et al.l . whose 
data was also based on UKIDSS-LAS, we see excellent agree- 
ment at the bright-end but a discrepancy in the faint-end 
slope. 

Fig. [S] explores the LF shape more closely by showing the 
1(7 error contours from our Schechter function fits in the 
M* — a plane, thereby illustrating the direction of the known 
degeneracy between M* and a. Recent values from the lit- 
erature, as indicated in Table [5] are shown as data points 
with error bars (colour coded on Fig. |B] according to fil- 
ter). The dashed lines on Fig. |B] are the ugriz contours 
produced when our sa mples are conservatively cut at the 
brighter limits used bv iMontero-Dorta fc Pradal . In general 
the optical data agree reasonably well with recent results 
from the t wo much larger SDSS studies of iBlanton et al.l 
(|2003f ) and lMontero-Dorta fc Pradal (|2009l ). and the conser- 
vatively cut sample and our standard sample 1 a M* —a con- 
tours overlap (except in the u band, which loses the largest 
fraction of gala^cies following the brighter apparent magni- 
tude cut). The volume of the common MGC-SDSS-UKIDSS 
region is much smaller, and the resulting uncertainty is sig- 
nificantly larger than these two previous SDSS results. As 
an aside, we note that the two SDSS results, while on the 
whole are consistent with our results, appear to be incon- 
sistent with each other at a high level of significance. After 
estimating the i?(Z)-uncertainty, it seems unlikely that this 
is the only cause of discrepanc y between the two SDSS re- 
sults (|Montero-Dorta fc Pradal do not use an evolution cor- 
rection, Blanton et al.l does). This suggests that significant 



unquantified systematics still remain and that the increase 
in statistical size from the MGC to the entire SDSS DR5 
dataset is not necessarily increasing the accuracy to which 
the measured LFs are known. 

The near-IR LF shapes, in comparison to previous stud- 
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Table 4. Derived Schechter function parameters in ugrizY JHK for the magnitude limits indicated within the redshift range 0.0033 < z < 
0.1. The errors shown for the Schechter parameters are, in order, due to the sample size, K(z)-correction and E(z)-corrcction uncertainties. 
The errors shown for luminosity density and ufi, statistics are due to sample size, and combined K+E correction uncertainties. These can 
be combined in quadrature to give the combined error, u and z results have been modified by —0.04 mag and 0.02 mag to compensate 
for the discrepancy between SDSS and AB magnitude systems. 



Sample (f>* (h^ Mpc~^) 



A/* (mag) 



UU15+U 
0016-0 
0011+0 
0008-0 
0011+0. 
0006-0 
0010+0 
0008-0 
0009+0. 

0011- 
0021+0 
0019-0 
0017+0. 
0016-0 
0013+0 

0012- 
0015+0 
0014-0 



UUU4+U. 
0004-0, 
0002+0. 

0002- 0. 
0001+0. 

0001- 0, 
0003+0. 

0003- 0. 
0002+0. 

0002- 0, 
0002+0. 
0002-0. 
0002+0. 
0002-0, 
0001+0. 
0001-0, 
0001+0. 
0001-0 



UU()3 
0003 
0001 
0001 
0002 
0002 
0002 
0002 
0001 
0001 
0003 
0003 
0002 
0002 
0001 
0001 
0000 
0000 



-18.2ltjj' 
-20.08io; 
-20.8llo 
-21.16lo; 
-21.46lo 

-22.20iQ 
-23.07lo 
-23.36i° 



U5+U. 

05- 0, 
05+0. 

06- 0. 
08+0. 

05- 0. 
07+0. 

06- 0. 
06+0. 

09- 0. 
12+0. 
11-0. 
11+0. 

10- 0. 
08+0. 
08-0. 
09+0. 
08-0. 



U4+U.03 
04-0.03 
06+0.02 
06-0.02 
03+0.03 
03-0.03 
01+0.01 
01+0.01 
01+0.01 
01-0.01 
01+0.01 
01-0.01 
01+0.01 
01-0.01 
00+0.00 

00- 0.00 
01+0.00 

01- 0.00 



-0.93+0 

3+0. 
s+0. 

-0.90tp 
-0.96i°; 



03+0. 
04-0, 
03+0. 
02-0, 
04+0. 

02- 0, 
03+0. 

03- 0, 
03+0. 

04- 0, 
08+0. 
07-0, 
07+0. 
07-0, 
05+0. 

05- 0, 
06+0. 
05-0, 



Ui+U.Ui 
02-0.01 
00+0.00 
00-0.00 
01+0.01 
00-0.01 
01+0.01 
00+0.01 
01+0.00 

00- 0.00 
01+0.01 

01- 0.01 
01+0.01 
01-0.01 
00+0.00 
00-0.00 
00+0.00 
00-0.00 



~} (xlO** h 
Lq Mpc-3) 



+0.24+0.15 
-0.25-0.14 
+0.33+0.17 
0.34-0.15 
+0.37+0.14 
0.40-0.14 
0.47+0.13 
0.56-0.12 
+0.68+0.12 
0.47-0.11 
+ 1.08+0.16 
0.93-0.15 
+0.82+0.13 
-0.72-0.12 
-1.11+0.05 
-0.99-0.05 
-1.49+0.13 
-1.44-0.13 



uU (xlO^'' h 
W Mpc-^) 
^ gg+0.25+0.15 

5.31 
6.35 
6.99 
6.89 
6.41 
4.95 
5.65 
3.48 



ufi, (corrected, 
same units) 



.5V+U.35 
.50-0.32 
.35+0.68 
.04-0.64 
.60+0.62 
.45-0.58 
.83+0.49 
.52-0.47 
16+0.37 
67-0.36 
89+0.42 
.24-0.40 
68+0.26 
38-0.25 
42+0.07 
.16-0.07 
85+0.07 
74-0.07 



u <20.84 
g <19.81 
r <18.76 
i <18.34 
z <18.07 
Y <17.38 
J <16.89 
H <16.12 
K <15.67 



0.0279Io 
0.0158l:o 
0.0124lg 
0.0120lg 



0.0109: 
0.0146^ 



+0. 


+0 
0. 

0.0155lo 
0.0149lo 
0.01561° 



-21.94'' 



-1.153 
-1.18j 
-1.18j 
-1.18j 
-1.06"* 



-0.99"^ 



1.91 
2.17 
2.29 
2.66 
3.07 
3.24 
3.17 
5.38 
6.98 



0.22-0.14 
+0.80+0.41 
0.61-0.38 
1.03+0.40 
0.93-0.38 
1.25+0.33 

1.03- 0.32 
+ 1.53+0.26 
-1.18-0.25 
+2.14+0.31 

1.66-0.30 
1.28+0.20 

1.04- 0.19 
1.16+0.05 
0.95-0.05 

+0.74+0.06 
0.65-0.06 



4.50: 
8.97; 

9.9i; 

10.3' 



9.71+1 

-1 
-1 
-1 
-1 

+0 



8.663 
6.53j 
6.89j 
4.01:" 



ies, show significant offsets from the MGC-SDSS-UKIDSS 
data. This presumably reflects the quality of the underly- 
ing imaging data with the near-IR data rapidly evolving, 
in resolution and depth, from the very shallow 2MASS data 
through to the less shallow UKIDSS LAS. Unfortunately our 
sample of a few thousand galaxies is insuffici ent in size to 
enab le a full bi- variate brightness an alysis (e.g.. [Driver et aP 
[200i). However. [Smith et al.1 (|2009l ') have quantified the sur- 
face brightness limitations of the UKIDSS LAS data in the 
JiT-band. As the source data for our study is the same as that 
used by ISmith et al.l (|2009l ) it is reasonable to adopt their 
Jf-band limit for completeness of Mk — Slogjg^ ~ —17.5 
mag. Using the mean colours as indicated in Table. [3] we 
re-derive the Y JHK LPs within the revised ranges with no 
significant change to the Schechter paramaters. This implies 
that whilst surface brightness selection bias is a concern at 
some level it is unlikely to be affecting our fits which are 
dominated by systems at the bright-end of the LF. 



^ Smith et al (2008) 
Q Eke et al (2005) 




Blanton et al (2003) 
Q Montero-Dotta et al (2009) 
Baldly et al (2005) 



g Koctianeck et al (2001 ) 
□ Bell et al (2003)', 



i Coleet al (2001) 
^ Jones et al (2006) 



-23 



M„-51ogh (mag) 



4.2 The Cosmic Spectral Energy Distributions 

We can calculate the luminosity density (^a) in each fil- 
ter from the integration of the Schechter luminosity func- 
tion, using the parameters listed in Table [4] and Equation [1] 
We convert these measurements (in units of KLq ;^Mpc~'') 
to energy per frequency interval {hWyipc~'^5Hz) which is 
more useful when comparing data with significantly differing 
filter widths and for comparison to model SEDs. To make 
this conversion we use: 

t( \ ^ ■ in-0-4(Aff7i ,-34-10) /^N 

where A is the effective wavelength, jx is the total luminos- 
ity density and Mq ^ is the absolute magnitude of the Sun 
in the specified filter. The constant value (34.10) derives 
from the definition of th e AB magnitude s cale (where 
3631Jy equates to mag; lOke fc Gunnlll983l ). Note that 
the filter width technically should appear twice but cancels, 
i.e. to derive the total flux through the filter one should 
multiply by the filter width (in Hz), however to make a 
useful comparison it is more logical to show energy per (5Hz 
which requires dividing by the filter width (in Hz). Our 
values for j and vf{v) are in the final two columns of Table 



Figure 6. ugrizYJHK Icr error contours for the LF fits from 
Fig. [5] with LF fits and a uncertainties from equivalent surveys. 
Dashed contours are from sampl es that have undergone th e more 
conservative cuts introduced by iMontero-Dorta &: Pradal . SDSS 
equivalent survey points have been transformed from the filter 
system at 2; = 0.1 to the filter system at 2 = 0. Note that un- 
certainties due to K-l-E corrections are not included in the error 
contours, and they purely show the uncertainty in the chi-squared 
best fitting. 



[l] and values for previous surveys are in the final columns 
of Tabled 

Figure [7] shows the position of our total luminosity density 
values compared with previous measurements. The step- 
function seen between previous optical and NIR surveys is 
no longer apparent in our data (except for a slight decre- 
ment in J) perhaps suggesting that cosmic variance may 
indeed have played a part in this discrepancy. In general our 
data points are consistent with what has been published 
before but provides a relatively smooth distribution within 
a single survey suggesting that constructing the CSED from 
within a single survey volume is critically important. 
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Figure 5. ugrizY JHK luminosity distributions and fitted Schechter functions, with comparison lines for Schechter parameters from 
equivalent surveys. The coloured lines show the best fit Schechter function for ugriz samples that have undergone the more conservative 
cuts introduce d by [M ontero -Dorta & ; Pradaj . Poissoncan unce rtainties are shown for ea c h bin . It should also be noted that the absolute 
magnitudes in iBaldrv et al .1 l|2005l Vl Blanton et al ] l|2003l ') and iMontero-Dorta fc Pradal ll2009l 'l use SDSS pass bands that have been red- 
shifted to z = 0.1, and therefore have been k-corrected (and evolved, where applicable) back to 2: = 0. The iNorberg et al. I ll2002l') and 
[Driver et al comparison lines in the g ban d use the similar b j and B filters, respectiv ely, which have be en tr ansformed to the q 

band using the assumption that B — V = 0.94 from lLiske et al" H200I), and filter conversions in lLiske et al. 1 1I2OO3I) and lBlanton fc Roweij 
1I2OO7) '). 



The CSED values are still dust uncorrected and in Fig. [8] 
we show the pre- (solid symbols) and post- (open sym- 
bo ls) attenuate d value s. We adopt the prescription laid out 
in [Driver et al.l (|2008l ) to correct our cosmic SED; this re- 
sults in corrections of x2.27, xl.69, xl.56, xl.47, xl.41, 
xl.35, xl.32, xl.22, and xl.l5 in u, g,r,i, z,Y, J, H, and 
K respectively (black, solid symbols). For comparison, we 
also show dust corrections ca lculated using the p rescrip- 
tion laid out in Section 3.3 of lCalzetti et al.l (|2000l ): using 
Es{B — V) = 0.16 (from the same paper), this results in cor- 
rections of X2.46, X2.02, xl.69, xl.51, xl.39, xl.3, xl.22, 
xl.l3, and xl.06 in u, g,r,i, z,Y, J, H, and K respectively 



(red, solid symbols). Overlaid are three expectations derived 
for us by Stephen Wilkins (priv. comm) from various cosmic 
SFH-l-IMF combinations using the PEGASE code. The blue 
curve is based on the cosmic star-f o rmati on history (CSFH) 
assembled by iHopk ms fc BeacomI (|2006l ) from direct mea- 
surements reported in the literature but with a Salpeter 
IMF flattened below 0.5Mq . The brown curve uses the same 
CSFH but adopts the best fitting IMF of I Wilkins etal] 
(200^ , which has a high-mass slope slightly shallower than 
the typical Salpeter value (i.e., —2.15 rather than —2.35). 
Both appear to be inconsistent with our data, with the data 
fitting the blue curve in the optical and then tending towards 
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Table 5. Schcchter parameters and luminosity densities for optical and NIR surveys in the literature 
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Montero-Dorta fc Prada (2009) 
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t adjusted to effective filter rest wavelength for object at z = 0.1 

X adjusted 



(and then propagated 
to h = 1. 



through to the calculation of j and uf,^ ) 



the brown curve in the NIR. An intermediate solution could 
perhaps be found which would fit. The purple curve also 
adopts a Salpeter I MF but with t he cos mic star-formation 
history derived by IWilkins et al.l (|2008l ) . This is based on 
constraints from the evolution of the total stellar mass his- 
tory and which generally predicts a l ower star-formation 
rate a t higher redshift than reported in iHopk ins fc BeacomI 
(j2006h . However, like the brown curve, the magenta curve 
fails to fit the data at shorter wavelengths. As we are in 
the process of assembling a much larger dataset with signif- 
icantly smaller uncertainties (particularly diminishing the 
effects of cosmic variance) we defer a detailed comparison of 
the CSED, CSFR and IMF. 

Here we conclude that our values are broadly con s istent with 
the range of r esults reported in IWilkins et al] (|2008h and 
IWilkins et al.l (|2008l ). and that refined measurements of the 
CSED should provide useful additional constraints on the 
CSFH and IMF. Such improvements should arise via the 
following measures: 

(i) A larger statistical sample 

(ii) Deeper NIR photometry 

(iii) Matched photometry /deblended solutions across all 
filters 

(iv) A full bi-variate brightness distribution to model the 
selection bias 

(v) More sophisticated modelling of the K{z) and E{z) 
corrections 

These improvements are currently in progress within the 
ongoing Galaxy And Mass Assembly Survey (GAMA). 

5 SUMMARY 

We have combined data from the MGG, SDSS and UKIDSS 
LAS to produce a master catalogue of --^ 10k objects with 



ugrizY J H K photometry and 100 per cent redshift com- 
pleteness. After careful consideration of the colour bias 
and restricting our sample to z < 0.1 we produce nine 
pseudo-flux limited samples comprising of between 3328 to 
1589 galaxies, depending on the filter. Using a modified 
SWML method and least-x^ fitting to a Schechter func- 
tion parametrisation, we recover luminosity distributions 
and functions in all nine filters. The resulting LFs in the op- 
tical typically produce luminosity functions consistent with 
previous publications from the SDSS collaboration albeit 
with slightly higher normalisations. In the NIR we recover 
significantly distinct Schechter funtion parameters which we 
attribute to the rapid improvement in NIR technologies and 
fields-of-view. Finally we derive the cosmic spectral energy 
distribution from u to K which is shown to be smooth. Af- 
ter correcting for dust attenuation we find that our observed 
CSED agrees within the errors with that expected and in 
due course these observations should provide useful new con- 
straints on both the cosmic star-formation history and the 
universality of the initial mass function. Work is currently 
underway to significantly improve upon the results reported 
here via the Galaxy And Mass Assembly (GAMA) survey. 
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Figure 7. Cosmic energy ou tput from 0.1 to 3 fim. The model line shows the BED of a 1 3.2Gyr Sa- t ype g al axy from the spectral 



Eke ct al. ( 20 oii".|Jones et al. I 200&l.ltluang 
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The survey has been supported through grants from 
the Particle Physics and Astronomy Research Coun- 
cil (UK) and the Australian Research Council (AUS). 
The data and data products are publicly available from 
|http: / /www.eso.org/'^jliske/mgc/] or on request from J. 
Liske or S.P. Driver. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the U.S. 



Department of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http:/ /www.sdss.org/ , 
The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partici- 
pating Institutions are the American Museum of Natural 
History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve 
University, University of Chicago, Drexel University, 
Fermilab, the Institute for Advanced Study, the Japan 
Participation Group, Johns Hopkins University, the Joint 
Institute for Nuclear Astrophysics, the Kavli Institute for 
Particle Astrophysics and Cosmology, the Korean Scientist 
Group, the Chinese Academy of Sciences (LAMOST), Los 
Alamos National Laboratory, the Max-Planck-Institute 
for Astronomy (MPIA), the Max-Planck-Institute for 



The ugrizYJHK luminosity distributions 13 





WTH08 SFH, WHTToe Base IMF 

HLlOB SFH, WHTTOa Base IMF 

HBOB SFH. WHTT08 Best IMF 

9 MGC-SDSS-UKIDSS LAS (observed) 

O MGC-SDSS-UKIDSS LAS (dust corrected) 

O MGC-SDSS-UKIDSS LAS (Calzetti correction) 

0.6 0,7 O.B 0.9 1 2 3 

A(^m) 



Figure 8. The cosmic spectral energy distribution as derived 
from tlie combined MGC-SDSS-UKIDSS LAS surveys com- 
pared to that expected from various cosmic star-formation his- 
tory/initial mass function combinations as indicated (Wilkins 
priv. comm.). Solid symbols are the raw empirical data un- 
corrected for dust attenuation an d open symbols are cor- 
rec ted for dust attenua tion following [Driver et al.l |2008| (black) , 
or ICalzetti et aD 120001 (red). Errors are only shown for the 
[Driver et al.l data for clarity, and are split into two components; 
the black errorbars show the uncertainty due to chi-squarcd fit- 
ting, the orange errorbars show the uncertainty due to k and e 
corrections. 



Astrophysics (MPA), New Mexico State University, Ohio 
State University, University of Pittsburgh, University of 
Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. 

The UKIDSS project is defined in lLawrence et al.|[2007l . 



Casah et al.l 


2007) and a 


Hewett et al. 


l2006l. The 



and science 

archi ve are described in Irwin et al (2008) and lHamblv et aP 
l2008l . 

STILTS is a set of command-line based catalogue 
matching tools, based upon the STIL library. It is currently 
being supported by its author (Mark Taylor), and can be 
downloaded from , http://www.star. bris.ac.uk/~mbt/stilts/; 
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APPENDIX A: PROBLEMS WITH CASU 
SOURCE EXTRACTION SOFTWARE 

The WFCAM Science Archive (WSA. iHamblv et al.llioo^ ') 
is the storage facility for post-pipeline, calibrated UKIDSS 
data. It provides users with access to fits images and CASU- 
generated object catalogues for all five UKIDSS surveys, 
along with copies of catalogues from a number of other sur- 
veys, including the SDSS and the MGC. The creators of 
the archive have also included cross matched data tables 
that link objects in the different catalogues. The LAS-MGC 
cross match table (lasSourceXmgcDetection) contains the 
ID numbers for all MGC objects within 10 arcsec of a LAS 
source, and the ID of that LAS source. Using this tool, we 
downloaded catalogues matching MGC- Bright B band lumi- 
nosities to their counterpart LAS Y , J, H and K band lumi- 
nosities. We required objects where both the MGC and LAS 
objects were definitely galaxies (MGC CLASS = 1 and LAS 
pGalaxy > 0.9; this criteria is only used here to guarantee a 
galaxy-only sample, and is not used where completeness is 
important), not in an MGC exclusion region {INEXR — 0), 
with good photometry (QU ALITY ^ 2) and with no major 
err ors in the LAS ob servation {ppErrBits < 256; follow- 
ing [Smit^^et^l][2Q0^, this criteria removes all objects that 
lie within dither offsets, all possible crosstalk artifacts, ob- 
jects with bad pixels and objects close to saturation). We 
also specified that the distance between the galaxy's cen- 
troid in the two surveys had to be no more than 2.5 arcsec 
apart (lasSourceXmgcDetection. distanceMins < 2.5/60). 
We define A'ukidss as the K band Petrosian magnitude 




UKIDSS 



Figure A2. Tiic fraction of galaxies deblended as a function of 
^UKlDSSi where X is Y, J, H or K. The red vertical lines are our 
band dependent apparent magnitude limits , the blue ver t ical li ne 
is the K band apparent magnitude limit of lSmith et al. ] ||2009|') . 



(laaSource.kPetroMag) taken from the WSA dataset. 
Unfortunately, a problem arises when the colour distribu- 
tion of the sample is plotted (see the left image of Figure 
IA7|| . There are a signific ant number o f very bright objects 
that are extremely red. ISmith et al.l ()2009h noted this is- 
sue and have traced this problem to a catastrophic fault 
in the CASU deblending algorithm, which is causing de- 
blended galaxies to become significantly brighter than their 
parent object, in some cases by several mag. This is high- 
lighted in the left hand image of Figure lAll while the unde- 
blended bright galaxies (black points) are clustered around 
the -Bmgc — -?fuKiDSs = 4 mag line, the deblended bright 
galaxies (blue crosses) are spread across a much wider colour 
range. 

Unfortunately, the pre-deblended data are not output by 
the pipeline (when the CASU source extractor breaks up 
brighter galaxies, it does not generate any parameters for 
the parent object; Irwin et al, in preparation), so this fault 
is not trivial to correct or quantify. It is possible to re- 
move the deblended galaxies from the sample by filter- 
ing upon an error bit designated as the deblender flag 
{lasSource.kppErrBitsk.QyS)mmQlO = 0) but this would 
leave us with a bia sed dataset. This option was adopted by 
ISmith et aLl(|2009l) . but Figure |A2] shows our reservation; al- 
though only 3.5 per cent of galaxies in the unlimited K band 
sample, 3.2 per cent in H, 3.1 per cent in J and 10.8 per 
cent in Y are deblended, they are not uniformly distributed 
across the apparent magnitude range. 

We initially attempted to overcome the problem by using 
the SDSS optical colours to predict the Y , J, H and K band 
fluxes for these problem galaxies. Unfortunately, the corre- 
lation was too noisy (Am ~ ±lmag), so we were forced to 
abandon the CASU DR3PLUS catalogue products provided 
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Figure Al. Comparisons between the MGC B-K values for K values taken from the UKIDSS Survey data (i^uKlDSs)! a-nd from 
Sextracted UKIDSS images (A'mgc)- The dotted green lines are the apparent magnitude cuts used for the K band sample, and the red 
lines are the colour median and colour outlier lines for each sample. 



through the archive and entirely re-derive all Y , J, H and 
K photometry from the reduced fits images. 

Al Reanalysis of UKIDSS data using Sextractor 

The WSA contains a tool (MultiGetlmage) for extracting 
1 arcminute x 1 arcminute fits images of 500 objects when 
given a list of their coordinates. The RA and Dec of the 
MGC-Bright galaxies were input into this program and 
images were extracted (as was a list of the galaxies where 
no image was available because of limited coverage). 

The SExtracto r object extraction program 

(| Bert in fc Arnout3 1 19961 ) was used to extract objects 
from these images (taking zeropoint, pixel size, gain, seeing 
and background levels from the UKIDSS fits file headers), 
and the fluxes and aperture centres for objects in each 
image catalogued. Both elliptical aperture Petrosian and 
Best fluxes were calculated for each object (the SExtractor 
Best flux generally uses the Kron flux, except in crowded 
regions where it uses a corrected isophotal flux). As Bmgc 
magnitudes were also calculated using Best apertures, 
we adopt Best magnitudes for our NIR magnitudes (we 
define Kmgc as the SExtracted-BEST magnitude and 
^^MGC, Petrosian as the SExtracted-Petrosiau magnitude). 



Distances between the position of each newly SExtracted 
object and the centre of that image were calculated (using 
the catalogued aperture centres and axis sizes taken from 
fits file headers). As all extracted images were centred on 
the position of an MGC-Bright galaxy, the object that was 
extracted closest to the centre of each frame is assumed 
to be the match to that MGC galaxy. The apparent 
magnitude of this object is calculated from its flux (F), 
exposure time (t), extinction (Ext), zeropoint (Zpt) and 
airmass (secxmcan; the former is taken from the SExtractor 
created catalogue, the rest are taken from UKIDSS fits file 
headers), using: 

F 

m = Zpt ~ 2.51ogiQ( j) - Ext.(secxmcan - 1) (Al) 

Again, we exclude all NIR objects that are > 2.5 arcsec from 
the MGC galaxy's centre. The right hand graph in Figure 
lAll show how the Bmgc-Js'mgc distribution compares with 
that from WSA dataset. The online archive contains 58 de- 
blended galaxies with Bmgc — ^ukidss > 6 mag. When 
re-extracted, only 2 of those galaxies have -Bmgc — Kmgc 
in that range (one is a distant galaxy that will be removed 
from our luminosity function when we impose redshift limits 
and the other is just a very red galaxy). Figure [K3\ contains 
K band images of two of the deblended MGC galaxies, the 
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Figure A3. The position of tlie Xmgc Kron aperture (grey 
ellipse) and the central position of the ii'uKlDSS galaxy (black 
cross) for the deblended galaxies MGC65412 and MGC05276 

position of their Kmgc and iiTuKiDSS apertures and the lu- 
minosity returned using each method. Figure [X4l shows how 
ifuKiDSS, Kmgc and T^MGCPetrosian photometry compares. 
We note that the dominant error in galaxy photometry is 
typically due to the flux measurement method. We find for 
an individual galaxy a typical uncertainty of ±0.04mag be- 
tween the elliptical aperture BEST and Petrosian methods, 
and therefore adopt this as representative of our standard 
luminosity error. 

A 1.1 Near-IR galaxies with no MGC match 

It is possible that there are galaxies in the UKIDSS sur- 
vey that have no match to the MGC, despite being inside 
its area of coverage. While some of these objects may be 
misclassified by the automated classification system, a num- 
ber of these objects could potentially be particularly red 
galaxies (e.g. heavily dust attenuated or exceptionally high 
redshift galaxies) that have slipped below the Bmgc ~ 20 
mag threshold and must be recovered for the calculation 
of the appropriate NIR luminosity function. We extracted 
fits files for all UKIDSS galaxies {pGal > 0.9) with X < 

Smog, threshold — (SmGG, threshold — ^Survey )median . Wc USed 

SExtractor to extract objects from these images (using the 
same parameters as in subsection I Alf) . and calculated their 
luminosity using the Best aperture. We exclude d all objects 
that w ere not covered by MGC CCDs. Following lLiske et al.l 
(|2003l ). we also excluded all objects with stellaricity > 0.98. 
The remaining objects that lie above our apparent mag- 
nitude limits (see subsection 12. 5p were checked by eye, to 
remove any that looked like they were noise (or were other- 
wise suspect) that had been misclassified. We found a total 
of 100 NIR galaxies without a MGC match brighter than 



our K band apparent magnitude limit, 76 in H, 73 in J and 
42 galaxies brighter than our Y band sample. However, as 
these galaxies are still below our B magnitude limit, they 
are not used in the calculation of our luminosity function. 

A 1.2 Colour outliers 

Whilst our new catalogues produced by SExtractor contain 
a smaller number of colour outliers, some still remain and 
can be seen at the extremes in Figure IA5|) . The median 
and standard deviation of B — X was calculated, and we 
define colour outliers to have a colour > Median + 3a or 
< Median — 3cr. We have manually reexamined these objects 
and, in 11 galaxies in the K band, 19 in H, 1 in J and 1 in 
Y, the program is incorrectly deblending an object. In these 
cases we modify the SExtractor deblending parameters to 
achieve a consistent deblending outcome. 
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Figure A4. Comparisons between the K band galaxy luminosities using magnitudes taken from the UKIDSS Survey data (^i'uKlDSs)i 
and from SExtracted UKIDSS images using circular Petrosian (^^'mgc Petroaian) ^'Hd Best (i^MGc) methods. As in Figure fAll blue 
crosses are galaxies that have been flagged as deblended in the UKIDSS survey data, and black dots are those that have not. 




Figure A5. -BmgCj^mgc ~ ^mgc a-nd ^MGCi^MGC ~ -^MGC graphs for X=Y, J, H and K band data. The vertical lines in the left 
hand graphs are the apparent magnitude cuts. All galaxies above these lines are excluded from our luminosity function samples. 



